Searching for thermal signatures of persistent currents in normal metal rings 
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We introduce a calorimetric approach to probe persistent currents in normal metal rings. The heat 
capacity of a large ensemble of silver rings is measured by nanocalorimetry under a varying magnetic 
field at different temperatures (60 mK, 100 mK and 150 mK). Periodic oscillations versus magnetic 
field are detected in the phase signal of the temperature oscillations, though not in the amplitude 
(both of them directly linked to the heat capacity). The period of these oscillations ($o/2, with 
<J>o = h/e the magnetic flux quantum) and their evolution with temperature are in agreement with 
theoretical predictions. In contrast, the amplitude of the corresponding heat capacity oscillations 
(several A?b) is two orders of magnitude larger than predicted by theory. 

PACS numbers: 73.23.-b, 73.23.Ra, 68.65.-k, 64.70.Nd 



INTRODUCTION 

At very low temperatures T, a small isolated metal 
ring carries a weak but non-dissipative current that is 
not destroyed even by the presence of disorder in the 
ring [1] . It follows from very general arguments that this 
persistent current / is periodic in the magnetic flux $ 
threading the ring [2 , 3 : 

7($,T)= yj/ m (T)sin(27rm— J, (1) 

where $0 = h/e is the magnetic flux quantum. The 
typical magnitude of the persistent current in a ring of 
circumference L is / ~ eE^h/h^ where E^h = hD/L 2 is 
the Thouless energy and D is the diffusion coefficient de- 
termined by the scattering of electrons from the disorder 
in the ring. For L ~ 1 /im, we find I ~ 1 nA. The pre- 
cise value of / and the direction in which this current is 
flowing are determined by the configuration of impurities 
and are thus random. 

Due to the weak magnitude of the persistent current, 
only the first two harmonics m = 1, 2 in Eq. ([!]) were 
observed experimentally [4HT0]. They were also in the 
focus of numerous theoretical studies (see Refs. [TTHT8] 
for a representative selection of theoretical works). The 
ensemble of the existing literature indicates that the sta- 
tistical properties of the first harmonics I\ can be un- 
derstood within the model of noninteracting electrons, 
whereas the second harmonics I2 is dominated by inter- 
action effects. In the experiments performed on single 
rings [HE] or on small ensembles of identical rings [71 [TO], 
both the first and the second harmonics of I have been 



detected. Hence, the measured period of current oscilla- 
tion with $ is $0. The mean values of I\ and 1^ their 
variances, and possibly the full probability distributions 
can be measured [18] . On the contrary, if large ensembles 
of rings (up to 10 T rings in Ref. [4) are used [HE], aver- 
aging over all the rings of the ensemble inevitably takes 
place and one only has access to the average current (/). 
Because (Ii) =0 [19 , the measured signal is mainly due 
to the second harmonics (I2) ^ 0. The observed current 
then oscillates with a period 3>o/2. 

The properties of the first harmonics I\ of the persis- 
tent current are nowadays relatively well understood. If 
early observations were somewhat contradictory [5] [7], 
more recent works have demonstrated an impressive 
agreement between experiment and the noninteracting 
electron theory [9] [10] [18]. The second harmonics 
on the contrary, still represents a challenge for both the- 
ory and experiment. Even though a consensus exists on 
the importance of electron-electron interactions to ex- 
plain its properties, the calculation assuming repulsive 
interactions [I2j [14] yields the magnitude of (I2) which 
is significantly smaller than the one measured in the ex- 
periments. In addition, the paramagnetic response of 
the rings at small magnetic field predicted by the theory, 
disagrees with observations [4j [8] . Diamagnetic response 
can be obtained in a theory assuming attractive interac- 
tions [13 . However, interactions which are sufficiently 
strong to reproduce the experimentally observed values 
of (I2} would induce a transition to the superconduct- 
ing state at temperatures which are too high to be com- 
patible with known properties of (some of the) metals 
used in the persistent current experiments: copper, gold, 
and silver. Indeed these metals do not exhibit super- 




CbOOOOOOOOOO 

ooooooooooo 
ooooooooooo 
ooooooooooo 
ooooooooooo 
ooooooooooo 
ooooooooooo 



conductivity even at the lowest temperatures accessible 
experimentally (down to 0.1 mK PS]). A possible solu- 
tion to this problem has been recently proposed by Bary- 
Soroker et at. [17 j: a tiny amount of magnetic impurities 
can destroy the superconductivity but has little impact 
on the persistent current. Two new parameters — the 
spin-scattering rate l/r s and the bare superconducting 
transition temperature T® of the material without mag- 
netic impurities — appear in the theory and allow for 
a reasonable explanation of the experiments reported in 
Refs. [4 and [7]. 

Motivated by the recent progress in the research on 
persistent currents in normal metal rings, we propose 
here a new way of detecting these currents. Our approach 
is radically different from the one employed in all previous 
experiments: if these anterior works relied on the mea- 
surement of the rings' magnetic moment (either using a 
more or less sophisticated version of a SQUID magne- 
tometer [4j |5j [71 [9], coupling the rings to a supercon- 
ducting microresonator |8 , or using an elegant microme- 
chanical detector [10 j), we propose to focus on the rings' 
heat capacity. Our idea stems from the basic principles 
of thermodynamics. On the one hand, the heat capacity 
at constant pressure is given by C p = —T(d 2 F/dT 2 ) pi 
where F is the thermodynamic free energy. On the other 
hand, the persistent current is / = —dF/d&. Therefore, 
we find 



dC 



dT 2 ' 



(2) 



From here on we omit the subscript 'p' of the heat capac- 
ity to lighten the notation. Equations ([!]) and Q imply 
that 



C($,T) = C(0,T) 

oo 



1 — cos 2iim 



, (3) 



where C m (T) = (T§Q/2imi)d 2 I m (T) /dT 2 . The heat ca- 
pacity C is therefore also a periodic function of the mag- 
netic flux 

Equation ([2| shows that the dependence of C on $ is 
intimately related to the dependence of I on T. This 
link between C and I was already exploited by Yang and 
Zhou to investigate the impact of spin-orbit coupling on 
both the persistent current and the heat capacity in the 
model of noninteracting electrons [20], as well as by (some 
of the) present authors to study the entrance of magnetic 
vortices into superconducting loops [21-23]. The exper- 
iments reported here are, however, much more involved 
than those performed in the superconducting state: not 
only we work at much lower temperatures (T ~ 100 mK 
instead of T ~ 1 K), but also the variations of the heat 
capacity to detect are much weaker (AC ~ fc# instead of 
AC~10 3 k B ). 




FIG. 1: (a) Scanning electron microscope (SEM) image of a 
single silver mesoscopic ring, the scale bar represents 200 nm. 
(b) SEM image of the array of rings. 



EXPERIMENTAL SETUP 

Our sample (Fig.JT]) is composed of N = 5 x 10 6 nonin- 
teracting silver rings (2r = 600 nm diameter, w ~ 140 nm 
arm width, d = 34 nm thickness, total mass m = 330 ng). 
The silver rings are deposited by e-gun evaporation un- 
der 10 -6 mb vacuum. Thanks to the high purity of the 
silver the phase coherence length is large, of about 10 
/im at the temperature of 100 mK, which is far above 
the diffusive ring circumference L [24] . The specific heat 
of silver at T = 0.1 K can be estimated to be of the order 
of 10 -6 J/gK, and hence the magnetic-field-independent 
part of the heat capacity of the rings is expected to be 
3.3 x 10 -13 J/K; representing approximately 10% of the 
total heat capacity. The rings have been patterned by 
electron beam lithography on the thermal sensor made 
of a suspended silicon membrane (size 4 mm x 4 mm, 
thickness ~ 5 /im). The transducers, a copper heater 
and a highly sensitive NbN thermometer are integrated 
on each side of the membrane (see Fig.|2|, the silver rings 
being located between the two elements [25] . 

A sketch of the ac calorimetry technique that we use is 
given in Fig. [2] and described in detail in previous articles 
[26H28] . It consists in applying an ac current through the 
heater; here at a frequency f e i ec = 30 Hz. This induces 
oscillations T ac (t) of the temperature of the suspended 
membrane. These oscillations are detected by the ther- 
mometer. For a specific experimentally determined oper- 
ating frequency, quasi-adiabatic conditions are fulfilled, 
allowing measurements of the specific heat. The typical 
signal obtained from the lock-in amplifier is composed of 
a modulus measuring a root mean square (RMS) voltage 
proportional to the temperature oscillation of the mem- 
brane (T^ MS oc ^ a ^ MS ) an d of a phase <p, both related to 
the heat capacity through the equations: 



Tac(t) = 

tan(A^) = 



ujC 
2ujC 



cos(2co>£ + (/?), 



(4) 
(5) 



where Po is the power dissipated in the heater, uj = 
2?r feiec the electrical excitation frequency, C the heat ca- 
pacity of the membrane and K the thermal conductance 
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FIG. 2: Schematic drawing of the experimental set-up of heat 
capacity measurement on a suspended silicon membrane. The 
internal oscillation signal of the lock-in ftl is used to pilot 
the voltage current convertor used to apply the current to 
the heater. The lock-in ft 2 is used to measure the resistance 
of the heater, and then estimate the power dissipated. The 
oscillation of temperature T ac is measured on the lock-in ftl 
at the second harmonic (2f). The scale bar represents 1 mm. 



between the membrane and the heat bath (for more de- 
tails see [ED ESI E3 ES]). The setup is cooled down to 
very low temperatures using a dilution fridge, equipped 
with a superconducting coil supplying a magnetic field 
H normal to the plane of the rings. As compared to pre- 
vious measurements [21-23 , the thermometry has been 
adapted to work at very low temperature and the arm 
width has been increased from 40 /im to 150 /am to al- 
low measurements at a higher frequency. Consequently, 
ST dc = T membrane - T cryost at decreased and we were able 
to obtain temperatures compatible with the purpose of 
our study (T ~ 70 mK on the membrane for a regulated 
temperature of 50 mK of the sample holder) [29] . 



EXPERIMENTAL RESULTS 

Figure [3] shows a scan of the heat capacity versus the 
temperature of the membrane. The heat capacity mea- 
sured at the lowest temperature is around 10 pJ/K, with 
the error not exceeding 0.1 pJ/K. It must be noticed 
that this value corresponds to the heat capacity of the 
whole sample (rings + membrane). In the inset of Fig.[3j 
we show a histogram of the noise measurement obtained 
from the upper panel of this figure, indicating a noise of 
~ 0.8 nV/vfe. 
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FIG. 3: (color online) Top panel: voltage at twice the ex- 
citation frequency versus time at T = 200 mK. Inset of the 
bottom panel: histogram of the noise extracted from the mea- 
surement presented in the top panel. Bottom panel: heat ca- 
pacity C of silver rings versus temperature in the absence of 
magnetic field. 



Our experimental method to detect persistent currents 
is based on scanning the heat capacity versus the applied 
magnetic field, at a constant temperature. During the 
measurement, two components of the temperature oscil- 
lation are recorded: the modulus and the phase. Both 
bring physical informations about the heat capacity vari- 
ations (see equations (4) and (5)). 

A typical scan is shown in Fig. [4] Starting at H = 
mT, we slowly increase the field to 50 mT by small steps 
(~ 1 mT). C is measured at each H for signal integration 
times of the order of 10 seconds. In order to improve 
the sensitivity and to detect the very weak oscillations, 
we realized a hundred of identical scans by carrying out 
cycles between H = —50 mT and H = 50 mT. These 
data are averaged during the signal processing described 
below. 

The measurements were performed at three different 
temperatures: T = 60, 100 and 150 mK. For each of 
them, the temperature of the sensor membrane oscillates 
with an amplitude of ~ 10% compared to the control 
temperature. 
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FIG. 4: Phase of the temperature oscillation (at the top) and 
modulus of the temperature oscillations, expressed in terms 
of heat capacity using Eq. (4) (at the bottom), of mesoscopic 
silver rings versus magnetic field measured at T — 150 mK. 

DATA PROCESSING 

In order to detect a periodicity in our data, a Fourier 
analysis of the phase of the temperature oscillation and of 
its amplitude, expressed in terms of heat capacity using 
Eq. (4), has been performed. As a first step, the low 
frequency trend observed in every scan has been removed 
with a third-order polynomial regression. As a second 
step, the power spectral densities have been calculated 
as follows. The autocorrelation function of the signal 
x(H) (where x is the modulus or the phase), 

1 f M 

r x (H) = lim — / x(H')x(H' - H)dH', (6) 

M-too ZM J _ M 

is first estimated for every scan. Averaging over a hun- 
dred of scans is then performed to reduce the noise. 
The power spectral density S x (v) is calculated as the 
Fourier transform of the averaged autocorrelation func- 
tion Y X {H). The power spectral densities (refered to PSD 
in the following) of both the phase of the temperature os- 
cillation and the modulus of the heat capacity (extracted 
from equations 4) for the three temperatures 60, 100 and 
150 mK are displayed in Fig. [5] . 

We first stress that the high peaks observed for all 
curves at a low frequency v ~ 0.05 mT" 1 are a remi- 
niscence of the trend which has not been completely re- 
moved by the polynomial regression. These peaks do not 
carry any useful signal. Unfortunately, they can mask 
any physically interesting signatures of persistent cur- 
rents that might be present at low frequencies and would 
correspond to the first harmonics C\ of the heat capacity. 

In addition to the low-frequency peak, our spectral 
analysis of the phase of the temperature oscillation (up- 
per panel of Fig. pi reveals peaks at v ~ 0.11 mT -1 
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FIG. 5: (Color online) Power spectral densities of the aver- 
age signal for phase (upper panel (a)) and in the lower panel 
(b), the PSD of the modulus of the heat capacity (extracted 
from the equation 4). In the upper panel the red and the 
blue curves have been shifted from the green for clarity, and 
in the medium panel only the red curve has been shifted by 
10~ 30 (J/K) 2 .mT. The vertical straight lines delimit the fre- 
quency range in which signatures of persistent currents are 
expected. 



(the period is 9 mT). This value corresponds to a half of 
flux quantum through a ring of 264 nm in radius. Our 
rings have the inner radius of 220 nm and the outer ra- 
dius of 360 nm. Thus, the observed spectral peak is in 
the range of frequencies in which signatures of the second 
harmonics of heat capacity oscillations, C2, are expected. 
However, as it can be seen from the lower panel of Fig. [5j 
no clear signature of spectral peaks at v c± 0.11 mT -1 is 
observed in the modulus of heat capacity, although the 
phase and the modulus of the oscillating heat capacity 
signal are related through Eq. The reason for ob- 

serving a peak in the phase signal but not in the modulus 
signal most probably resides in a lower sensitivity on the 
modulus signal. At 100 mK, the sensitivity obtained for 
the modulus is - lO" 14 J/K or - lO" 20 J/K (~ 100 k B ) 
per ring (a sensitivity of AC/C ^ 5x 10 -4 ) as compared 
to the sensitivity obtained for the phase estimated to be 
of the order of A(p/(p ~3x 10 -5 . 

In order to compare our measurements with theoretical 
predictions, we estimated the average amplitude of the 
second harmonics of heat capacity oscillations, | (C2 (T)) | , 
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TABLE I: The amplitude of the second harmonics of heat capacity oscillations | (C2) | is extracted from the curves of Fig. [5] 
(upper panel). First, the area A of the spectral peak is estimated (second column). For this purpose, we fit the curves of the 
upper panel of Fig. [5] with smooth polynomials and evaluate the integral of the difference between the data and the fit taken 
between the vertical lines. Next, the associated heat capacity is calculated using Eq. ^ with <p = \/2A (column 5). Finally, 
the heat capacity is divided by the number of rings iV = 5x 10 6 to obtain the signal per ring (last column). 



from the area below the curves of Fig. [5] (see Table [I]). 



COMPARISON WITH THEORY 

We compare our results with the predictions of two 
theoretical models: the model of Ambegaokar and Eck- 
ern (AE) [12] and the model of Bary-Soroker, Entin- 
Wohlman and Imry (BEI) [17]. Both models rely on 
taking into account interactions between electrons in the 
rings, but the interactions are repulsive in the AE model 
and attractive in the BEI model. In order to suppress 
the superconductivity that strong attractive interactions 
may induce, the latter model includes scattering of elec- 
trons by magnetic impurities. 

Even if the applicability of the AE model to realistic 
experiments is made questionable by the disagreement of 
the sign of the current that it predicts with measurements 
[8], it is still of great interest to describe experimental 
data [10 . Using Eq. (18) of Ref . [12] and our Eq. we 
readily obtain: 



<C 2 AE (T)> 
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FIG. 6: (Color online) Average amplitude of the second har- 
monic C2 of the heat capacity oscillations with the magnetic 
field. Points show the experimental results from Table [I] with 
the scale on the left. Lines show theoretical predictions: Eq. 
([7]) (dashed blue line) and Eq. (|8| (solid red line), scale on 
the right. 



where N(0) is the electronic density of states at Fermi 
energy and V is the mean value of the attractive in- 
teraction potential [12]. Using the data from Ref. [8], 
where silver rings similar to ours were studied, we esti- 
mate N(0)V ~ 0.37 and E Th /k B ~ 40 mK. The resulting 
dependence of (C^ E ) on temperature is shown in Fig. pi 
by a dashed line, with the scale given on the right. The 
theoretical Eq. ^ reproduces the trend of the temper- 
ature dependence of our data that has a maximum at 
T ~ 3E Th /k B « 100 mK. But the values of |(C^ E )| are 
2 orders of magnitude smaller than the data. 

In the BEI model, assuming attractive interactions in 
combination with scattering on magnetic impurities [17] , 
we find: 



<C 2 BEI (T)> 



d 2 



77. * 



dx sin(27rx) 



ty'[Fn(x, T, T c °, T s )] 



9T 2 1 ln(T/T°) + *[F n (x, T, T°, t s )} - *(l/2) . 



(8) 



where Sk(x) and ^ f (x) are the digamma function and 
its derivative, respectively, r s is the spin-scattering time, 
is the bare superconducting transition temperature in 
the absence of magnetic impurities (i.e. at r s — » 00), and 
F n (x,T,T°, Ts ) = (l + \n\)/2 + (l/2nT s + nx 2 E Th /k B )/T. 
The minus sign in Eq. ([8| reflects the diamagnetic na- 
ture of the persistent current in this model, in contrast 
to the paramagnetic current leading to Eq. (|7|). It fol- 
lows from the comparison of the theory with previous 
measurements of persistent currents [4j [7] that it is rea- 
sonable to choose T c ° = O.l^Th/te HZ]. At the same 
time, s = l/7rT^r s should exceed 0.87 to suppress the 
superconductivity at all temperatures. We set s = 1 and 
show the absolute value of the resulting heat capacity 
I (C 2 BEI ) I in Fig. H by a solid line. 
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DISCUSSION 

Interestingly, the two theoretical approaches consid- 
ered above yield very similar results both for the mag- 
nitude of heat capacity oscillations ( | (C2 ) | ~ 10 -2 /cb) 
and for its temperature dependence (maximum values of 
|(C2)| reached at T « 100 mK for experimental parame- 
ters). Whereas the latter temperature dependence is in 
agreement with our measurements, the predicted magni- 
tude of the heat capacity oscillations is far too small to 
explain the observed | (C2) | °° &b- This discrepancy may 
result from the insufficiencies of the theoretical models, 
as well as from the uncontrolled errors in the estimation 
of the absolute values of heat capacities in our experi- 
ment. 

On the theory side, the two models considered here 
were previously shown to be compatible with experiments 
(see Refs. [12] and [H]). They provided reasonable re- 
sults for the magnitude and the temperature dependence 
of persistent currents measured in Refs. [4] and [7]. It 
should be noted, however, that, on the one hand, the 
heat capacity is proportional to the second derivative of 
the persistent current / with respect to temperature [see 
Eq. ([2])] and hence, it is sensitive to fine details of the 
temperature dependence of / that might not be captured 
by the theory. Equation ([8|, for example, changes sign 
for T « 30 mK, leading to a peculiar low-temperature 
behavior of the heat capacity in Fig. [6j This is due to 
the change of the shape of {^(T)) curve from convex 
to concave. Also the surface states (evanescent states), 
when included in the theoretical model, may modify its 
predictions significantly [30]. In addition, it was noticed 
by several authors that other phenomena (such as, e.g., 
the fluctuations of electron spin density [31] or the ambi- 
ent electromagnetic field [32]) may induce magnetic flux- 
periodic currents in mesoscopic rings. These currents can 
be comparable or even larger than the persistent currents. 
Given the large signals measured in our experiments and 
taking into account that our calorimetric technique may 
be particularly sensitive to currents that dissipate heat, 
we believe that it is likely that these phenomena may be 
important in our setup. Because Eq. ([2| does not hold 
for these non-equilibrium processes, it remains to be seen 
if and how they could be included into the theoretical 
model. 

The disagreement between theory and experiment seen 
in Fig. [6] might also stem from the difficulties in obtaining 
quantitatively correct values of heat capacity, intrinsic to 
the extreme difficulty of the measurements: measuring 
very low thermal signals at extremely low temperatures. 
These difficulties are obvious already from the compar- 
ison of the heat capacities extracted from the absolute 
value and the phase of the voltage signal measured in 
our experiment (compare the two panels of Fig. |5|. 



CONCLUSION 

We introduced a calorimetric approach to the study 
of persistent currents in mesoscopic rings made of nor- 
mal metals. The approach relies on the measurement 
of periodic variations of heat capacity of a large ensem- 
ble of rings with magnetic field. Under the experimental 
conditions reported here, the approach was at the limit 
of its sensitivity, the signal being strongly masked by 
noise. Despite this, we estimated the amplitude of heat 
capacity oscillations to be of the order of several &b per 
ring at ~ 100 mK. The amplitude is two orders of mag- 
nitude larger than expected from the existing theories 
which, however, correctly predict the range of tempera- 
tures where the heat capacity signal is maximum. Fur- 
ther experiments and theoretical investigations are nec- 
essary to elucidate the sources of this discrepancy. 

Both experiment and theory suggest that, in contrast 
to the persistent current /, the heat capacity C of an en- 
semble of mesoscopic rings is not a monotonic function of 
temperature. In particular, the average value of its sec- 
ond harmonic C2 vanishes for T — >• 0. More experiments 
having better sensitivity will be necessary to evidence the 
real position of this maximum. 
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